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I. INTRODUCTION 

The problem of dissipative dynamics in quantum system is a key issue in physics and chemistry 
today. The system-reservoir model describing the evolution of the quantum system coupled to 
a much larger system regarded as a reservoir has been the standard paradigm for dissipation 
in classical and quantum systems for many years [1-3]. A popular variant of this model is the 
spin-boson model, a magnetic dipole coupled to a Boson field; another is a two-level atom in 
contact with a continuum of radiation field modes. These models describe a large variety of 
physical situations, such as, spontaneous emission, polaron formations, exiton motion, macroscopic 
quantum tunneling etc., in atomic physics, solid state physics and quantum optics. 

In general, two very distinct situations emerge depending on the strength of the coupling 
constant between the system and the reservoir. In the weak coupling case, the behavior of the 
former is only slightly affected by the reservoir which essentially behaves as a free field. In 
the strong coupling cases (as in polaron theories) the polarization of the reservoir field by the 
system can not be ignored. The second important approximation that is almost always made 
is that the correlation time of the reservoir must be very short enough (Markov approximation) 
for the interaction between the system and the bath to be small (weak coupling approximation). 
Although several generalizations of the theory [3] which go beyond these approximation schemes 
are now available which describe the various interesting physical situations in condensed matter 
and quantum optical physics, one essential step is the assumption of an equilibrium distribution 
of the reservoir modes. Very little attention has been paid to the problem where the reservoir in 
contact with the system is itself not in equilibrium. This nonstationarity of the bath is known to 
affect the kinetics of the system [4-6] leading to nonexponential decay in contrast to exponential 
decay in equilibrium activation rate theories. Thus the relaxation of the nonequilibrium modes 
may influence the dissipation of the system in question in a nontrivial way. The present paper 
addresses a related issue pertaining to quantum optical situations. 

We extend the quantum theory of dissipation of a harmonic oscillator coupled to a bath where 



the bath in question is in a nonequihbrum condition. The nonequihbrium bath is effectively 
reahzed in terms of a semi-infinite dimensional broad-band reservoir which itself kept in contact 
with a thermal reservoir. We make use of a systematic time scale separation to construct the 
appropriate Langevin dynamics of the system mode. The fluctuation-dissipation and the Einstein 
relations have been suitably generalized. A detailed study of two model cases as immediate 
applications has been carried out. We show that the Wigner-Weisskopf decay rate constant of the 
system mode is dynamically modified when the system is coupled to a set of relaxing modes. We 
also calculate the transient noise-spectra of the cavity mode with a positive gain. The result is 
remarkably different from the steady-state spectra. 

The outline of the paper is as follows: In section II we first generalize the quantum theory of 
dissipation of a harmonic oscillator for a bath which is not in thermal equilibrium, followed by a 
derivation of fluctuation-dissipation theorem in the next section III. Section IV is devoted to the 
discussion of the Wigner-Weisskopf decay of the system. In section V we calculate the transient 
noise spectra of the system mode with a positive gain. The paper is concluded in section VI. 



II. THE MODEL AND THE EQUATIONS OF MOTION 

To start with we consider a model consisting of a harmonic oscillator (the system) coupled to a 
set of relaxing modes considered as a semi-infinite dimensional system which effectively constitutes 
a nonequihbrium reservoir. This in turn is in contact with a thermally equilibrated reservoir. Both 
the reservoirs are composed of two sets of harmonic oscillators characterized by the frequency sets 
{ujj} and {^j} for the equilibrium and nonequihbrium bath, respectively. The system- reservoir 
combination develops in time under the influence of the total Hamiltonian 
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+ nj:9,iC,a^ + Cla)+hY.T.<^Ub]C, + b,Cl) . (1) 

The first term on the right-hand side describes the system mode with characteristic frequency 
Uq . The second and the third term represent the thermal and the nonequihbrium hnear modes. 
The next two terms represent the couphng of the nonequihbrium bath with the system mode and 
the thermal bath where the coupling constants are g^ and aj^, respectively. In writing down the 
Hamiltonian we have made use of the rotating wave approximation. 

The Heisenberg equations of motion for the system and the reservoir operators at any given 
time is given by 

d(t) = -iuJoa{t) -i^g^C^{t) , (2) 

hj{t) = -iujjbj{t) - i J2 (^j^^C^^{t) , (3) 

C^it) = -in^C^{t) - ig^a{t) - I J2 ajM^) ■ (4) 

j 

Making use of the formal integral of the Eq.(3) for bj{t), 

in Eq.(4) we obtain 
C,{t) = -iQ^C.it) - ig^ait) - I Y. a,^e-'^^^'-'%{to) -^.11 «jm«.- f dt'C,{t')e-'^^^'-''\ (5) 
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Taking into consideration [1] that the interference time of I]j '^j>'^>^~*'^^*^*~* ■* ^^ niuch smaller 
than the time over which the significant phase and amplitude modulation of the linear modes 
C^{t) take place, the last term in Eq.(5) can be identified as a relaxation term in the usual way 
with damping constant 

7^, = 71 a^^in,) a.,(fi.) Din,), (6) 



where D[fii,) represents the density of states of the equihbrium modes evaluated at Q,y. Thus one 
can write down the Langevin equation of motion for the nonequihbrium mode C^ as follows; 

C,it) = -in,C,{t) - ig,a{t) - E7^.a(t) + f,{t). (7) 

Here the last term f^{t) represents the usual noise operator arising out of the coupling of the 
nonequihbrium modes with the thermal bath modes as given by 

/,(t) = -^^a,,e— ^■(*-*«)6,(to) , (8) 

3 

with reservoir average of /^(t) is zero, i.e., 

{f,{t))B = , (9) 

where by the average {0(t))^ of an operator 0(t) we mean (0(t))^ = Tr{0{t)pB}- Here pb 
denotes the initial density operator for the thermal bath {bj} and is a multimode extention of the 
usual thermal operator. This is given by 

where T is the equilibrium temperature. Note that in defining the average we assumed the initial 
factorization of the total density operator into subsystem densities for the system, thermal bath 
and the nonequihbrium bath. 

Introducing the slowly varying operator 

C,{t) = C^(t)e*^-(*-*«) (10) 

the Eq.(7) reduces to the following form; 

V 

where 

F^{t) = ^(t)e^^-(*-*») . (12) 



The relevant properties of the noise operator -F)j(t) can be summarized as; 

{FMb = > (13) 

and 

{Fl{t)F^it'))^ = %M^,)6{t - t')5,, . (14) 

The last relation follows from 

(6]„(to)&n(to))iJ = N{uJn)5rr,n , (15) 

where N{oOn) is the thermal average of the number operator of the equilibrium bath and is given 
by Niojn) = / Jn s j - . Also note that 5^y takes care secular approximation. Eq.(14) also implies 
a purely Ohmic frequency-independent dissipation of the nonequilibrium modes. 

Taking into consideration of the standard fluctuation-dissipation relation for the thermal bath 
in terms of Eq.(14), we arrive at the following Langevin equation for the nonequilibrium bath 
modes; 

C,{t) = -tn,C,it) - tg,a{t) - ll,C,{t) + f,{t) . (16) 

Eq.(16) constitutes an important result of this section which takes into account of the relaxation 
of the intermediate oscillator modes due to their coupling to the standard thermal bath whose 
fluctuations are described by /^(t). It is important to note that the above consideration is based 
on the rotating wave approximations (RWA) which is frequently used when considering coupling to 
a heat bath in quantum optics. In the present case as we are to see in the subsequent sections that 
there are quite subtle interaction effects. The question whether these effects would survive a more 
complete treatment may naturally arise. We mention two pertinent points at this stage. First, 
if one retains the nonrotating couplings in the Hamiltonian and carry out the same perturbative 
procedure one arrives at an equation of motion for C^{t) [instead of Eq.(ll) ] which additionally 
contains secular oscillating terms of the form (7tg«(f^M+f^i')(*-*o) ^ KWA amounts to neglecting these 
contributions which may be important only at very high coupling. Secondly, C^ modes executes 
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a slow relaxation dynamics on the time scale ~ 1/7hu ^^ compared to the time scale of correlation 
of thermal noise. The time scale of secular oscillations being short, those can be safely averaged 
out from the relevant dynamics. 

Another important point regarding RWA in the context of the present linear coupling model 
Hamiltonian (1) needs to be considered. Ford, O'connell and Lewis [13] have demonstrated that 
independent oscillator model within RWA (a variant of LC model) is seriously flawed since the 
Hamiltonian becomes imaginary when the bath is not passive ( i.e., there exists an associated 
spectrum of eigenvalues ranging upto — oo ). The problem essentially lies at the specific frequency 
dependence of the coupling constant in its denominator (for example, as shown in Ref. [13], 
appropriately transformed Qj in Eq. (1) is proportional to J^ ). However, a standard procedure 
in quantum optics is to replace the summation over modes by an integral over their density D{ujj), 
(in free space equal to {-^)), i.e., Yjj{9jY -^ I-^ dujD{uj)g'^{uj) . One thus gets rid of the 
unwanted frequency dependence in the denominator of the coupling constant in the calculations 
(see, for example, calculation of Wigner-Weisskopf decay rate. Lamb et al in Ref.[l] ). Since we 
have followed the same approach, RWA does not pose any problem in the present analysis. The 
problem however, remains for a strictly discrete spectrum. 

We have presented above an extension of the quantum theory of damping from Langevin point 
of view within a traditional system-reservoir linear coupling scheme. Essentially the model consists 
of replacing the reservoir by damping terms in the Heisenberg equations of motion for dissipation- 
free system and adding fluctuating forces as driving terms which add fluctuations to the system. 
The operator forces are such that (i) the system has the correct statistical properties to agree 
in the classical limit and (ii) they maintain the commutation relations for Boson operators to 
ensure that uncertainty principle is not violated. These considerations have been fully taken care 
of in our analysis with appropriate elaboration in the following sections. The spiritual root of 
the quantum statistical approach to damping lies in the fluctuation-dissipation relation, which 
illustrates a dynamical balance of inward flow of energy due to fluctuations from the reservoir into 
the system and the outward flow of energy from the system to the reservoir due to the dissipation 
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of the system mode. We address this specific issue in the next section. 

III. FLUCTUATION-DISSIPATION RELATION FOR NONEQUILIBRIUM BATH 

To explore the influence of an initial excitation of the semi-infinite dimensional intermediate 
reservoir modes and its relaxation, we now consider the evolution of these linear modes C^ in 
terms of Eq.(16). The physical situation that has been addressed is the following; 

We consider that at t = to the excitation is switched on and the bath modes (C^, C^) are 
thrown into a nonstationary state such that they behave as a nonequilibrium reservoir undergoing 
relaxation. We follow the stochastic dynamics of the system mode and the relaxation of the 
nonequilibrium reservoir modes after t > t^. We assume that the effect of back reaction of the 
system mode on the reservoir modes is small enough to be neglected. Eq.(16) allows a formal 
solution of the following form 

J to 

The first term on the right hand side in the absence of the coupling of the system mode 
represents the ( long time ) stationary stochastic solution of the form 

C^(t) = c'^^-i[^^.{t-tn)+r^] ^ (18) 

where the amplitude C* (operators) and phases (j)^ (c-numbers) are assumed to be randomly dis- 
tributed. The random nature of C^(t) may be understood in the following way : Let us first note 
that in the absence of coupling g^ Eq.(16) allows the solution 



J to 



In the steady state we neglect the first term which decays rapidly. The second term is a 
randomly fluctuating term (which is the most important term in almost any Langevin analysis) 
due to f^(t). Substituting Eq.(8) in the above equation we obtain 



, +\ =p-^^M(t-to) V('-7n....^ p-7^p(t-to) hJfr.\ f 



Clit) = e"^^^(*"*°)^(-mj^) e-^^^(*-*°) 6j(to) / dt' e(^^^+^f^M+7^M )(*'-*«) ^ 



where s signifies the steady state. The above solution implies that C^(t) is essentially a super- 
position of unknown (since the initial condition for the infinite number of thermal bath oscillator 
hjito) are assumed to be completely uncertain) amplitudes (operators) and phases (c- numbers) 
and may be written compactly in the form of Eq.(18). When written in the form (18) we obtain 
an instantaneous realization of the random distribution of C^ and 0fj. Thus Eq.(17) represents an 
instantaneous solution of Eq.(16). To check the consistency of the solution (17) we note the follow- 
ing points; (i) To recover initial deterministic amplitude of C^{t) we need an ensemble average of 
Eq.(17) at t = to so that ^^^(to) = 0. Also we are to see (subsequent to Eq.(24)) that the first term 
in Eq.(17) is responsible for the usual fluctuation-dissipation relation when the nonequilibrium 
modes get equilibrated at t = oo . CHt) in C^(t) is thus dictated by the condition of stationarity. 
(iii) It can be easily seen that the random nature of the operator forces in CHt) is responsible for 
maintaining the Boson commutation relation for C^{t) which further ensures that the uncertainty 
principle is not violated. 

The second term on the right hand side in Eq.(17) carries the information of relaxation of the 
C^ modes due to their coupling to the thermal bath and is a typical memory term. The third 
term on the other hand represents the effect of coupling of the system mode to the nonequilibrium 
reservoir. 

Substitution of the solution(17) in Eq.(2) yields the equation of motion for the slowly varying 
system operator A{t)] 



i(t) = -y2gl [ rft'e*("«-^'')(*-*')e-^MM(*-*')^(t') + z{t), (19) 



^ "'to 



where 



A{t) = a(t)e-*^"(*-*«) (20) 



and 



Z{t) = -tY.9t^[C!.{t) + C^(to)e(-*^''-^MM)(*-*o)]e-o(t-to), (2I) 






Eq.(19) is a non-Markovian equation where the memory effects arise out of the two sources; 
the first one being the system operator concerned term A{t') which depends on earher time t'; the 
other one e"'^^'''''*^* •* is due to the relaxation of the nonequihbrium bath modes which arises because 
of their couphng to the thermal bath. Z{t) represents the noise operator for the nonequihbrium 
bath modes. 

In the weak coupling approximation scheme the first term in equation(19) can be simplified to 
the following form; 

A{t) fdnp{n)g\n) /°° rfre^(^"-^)"e-^MM-, 





where the summation over the bath modes is replaced by integration and p(f2) represents (an 
a priory known) density of nonequihbrium bath modes. Assuming weak dependence of 7fj on 
the modes one can reduce the last expression (after performing integration over Q) to obtain the 
following Langevin equation for the system operator; 

A{t) = -T A{t) + Z{t), (22) 
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where 



T = 'Kg\u^)p{u^) (23) 

can be identified as a dissipation constant of tlie system mode due to tlie fluctuations of tlie 
nonequilibrium reservoir in tlie limit when 7^^ is vanishingly small. The explicit dependence of 
the decay constant of the system mode on 7 will be revealed clearly ( as a more general case ) in 
the next section within the scope of Wigner-Weisskopf approximation. 

The nonequilibrium generalization of the fluctuation-dissipation relation is now immediately 
apparent. From the expression for Z{t) (Eq.(21)) one finds that 

We denote the average photon number of the nonequilibrium bath by 

niQ„to) = {Cl{to)CM)^^ , (25) 

where to signifies the dependence of average photon number of the nonequilibrium bath on its 
initial state of preparation. Also the steady state average photon number is given by 



By {0{t))iyji we mean {0(t))iyji = Tr{0{t)pc} where pc indicates the initial thermalized density 
operator for the intermediate oscillator {C}-modes, and is given by 

Pc = Y[exp{ ^^^ '" )[!- exp{^)] . 
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As usual, the initial factorization of densities of {b} and {C} modes is assumed. 

After replacing the summation by integration and 7^ by the average 7 in Eq.(24) we obtain 

Rewriting rn{ujo,to) in terms of a deviation from its steady state value rn(cjo) as 



Tn{uJo, to) = D{to) - Tn{uJo) , 
we identify a time-dependent diffusion coefficient D{t) in the last equation as 

Dit) = Tniuo) + [Dito) - rn(cuo)] g-^^^*-*") . 
We thus obtain 

(Zt(t)Z(t'))NR = {Tniuo) + [D{to) - rn(cuo)] e-^^^*"*")} (5(t - t') . (26) 

Eq.(26) is the desired nonequilibrium generalization of the ffuctuation-dissipation relationship. 
This relates instantaneous ffuctuations of the nonequilibrium bath (which itself is undergoing re- 
laxation at a rate 7 due to its coupling with the thermal bath) to the dissipation of the energy 
of the system mode through F. The nonequilibrium nature of the bath is implicit in the initial 
preparation which creates an initial diffusion coefficient -D(to) and also in the exponentially de- 
caying term. In the long time limit one recovers the usual ffuctuation-dissipation relation for the 
thermal bath at equilibrium. 
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We now express Eq. (26) in terms of energy density of fluctuations of the nonequilibrium modes. 
The energy density which is proportional to the power spectrum centered around uo is given by 



u{^,t) = if/^ dT{Z\t)Z{t + r))e 



i(f2— tJo)r 



= -(]n(f]) + e-2''(*-*")[M(fi,to)--fin((])] . 

It is important to note that t is the slow time variable which is well separated from the time scale 
of thermal noise. The fluctuations of the noise operator Z{t) is now explicitly determined by the 
nonequilibrium state of the intermediate oscillator bath through its energy density -^(^2, t) at each 
instant of time t. In other words the instantaneous nonequilibrium energy density distribution of 
the fluctuating modes is related to the friction coefficient of these modes on the system degree of 
freedom through a dynamic equilibrium. One can immediately recover the classical version of the 
above equation in the high temperature limit ( where niVt) = — ^ — ~ 4r ) ^o obtain 



u{n,t) = -KT + e-2^(*-*«) \u{Q,to) - ^KT 



2 



2 



The above equation was derived earlier [6] in the context of classical relaxation kinetics of complex 
nonlinear systems. This reduction to classical version of the nonequilibrium fluctuation-dissipation 
relation serves as a consistency check for its quantum generalization (26) which is more relevant 
in quantum optical issues that we address in this paper. 

Another point should be emphasized regarding the fluctuation-dissipation relation (26). The 
very notion of a nonequilibrium bath apparently suggests that frequency distribution function of 
the modes of this bath should be a function of time and in principle, one should seek for it as 
a self-consistent solution from a quantum kinetic analysis. We have followed here an alternative 
route. Since we are working with Heisenberg operator equations of motion, the knowledge of initial 
total state density which is factorizable in subsystem densities (for the system, thermal bath and 
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nonequilibrium or intermediate bath ) at t = 0, is sufficient to describe the complete dynamics in 
terms of average values and correlation functions. We are concerned here (and also in classical 
theory [6]) with a priori given frequency distribution functions for equilibrium {bj} bath {D{uj)) 
and nonequilibrium {C^} bath (p(fi)) which are independent of time (nor we have introduced any 
non-thermal temperature in describing the nonequilibrium bath). The essential content of the 
nonequilibrium nature of the bath rests on a time-dependent energy density fluctuation distribu- 
tion function u{Q,t) as described above (varying over a slower time scale compared to the time 
scale of thermal noise ) which is a derived quantity rather than a self-consistently obtained func- 
tion which might be obtainable from a quantum kinetic analysis. The effect of initial excitation is 
to create an energy density function u{Q, 0) which differs from its equilibrium counterpart. This 
departure sets in a nonequilibrium situation. All these considerations also apply to the classical 
version of nonequilibrium fluctuation-dissipation relation [6] . 

Before bringing an end to this section some pertinent points are to be noted. First, we need 
to stress that in the derivation of the relation (26) we assume that Z{t) is effectively stationary 
on the fast correlation time scale of the thermal bath. Second, the theory as developed above is 
based on the consideration of quantum optical situations in mind. It is well known [2,11] that, in 
general, the noise from the equilibrium bath at low temperatures is very different from a simple 
white noise and concerns expressions which contain integral from distribution function of the 
bath over all frequencies. Most often such situations are encountered in condensed matter and 
in chemical physics of complex systems. However in the problems of quantum optics where the 
harmonic oscillator bath serves as a standard paradigm for optical fields one routinely uses a broad 
band white noise spectrum such that X)^ fi'u^M ^^ slowly varying and the summand in Eq. (24) is so 
strongly peaked at Q^ = Uq, that we may convert the sum into an integral and remove the slowly 
varying factors to obtain the result (26) [ p. 422 of Louisell in Ref. 1 ]. Similar consideration leads 
us to Einstein's spontaneous emission coefficient or Wigner-Weisskopf decay rate which contains 
single frequency t^o, the characteristic frequency of the system mode. Thus although, in principle, 
it may be possible to consider a colored noise spectrum and the resulting frequency dependence 
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of rate constants we restrict ourselves to former situations of broad band reservoir which itself is 
undergoing a relaxation on the time scale of I/7. An important content of the present work is to 
explore the effect of this secondary relaxation on the primary kinetics of the system mode. The 
nonequilibrium generalization of the fluctuation dissipation relation as discussed in this section 
serves as a basis of this exploration in following two quantum optical cases. 

IV. DECAY OF THE SYSTEM MODE ; WIGNER-WEISSKOPF 

APPROXIMATION 

We now obtain the solution of Heisenberg-Langevin integro-differential equation of motion 
[Eq.(19) ] for the system mode coupled to nonequilibrium bath under Wigner-Weisskopf approxi- 
mation. If we take the Laplace transformation of Eq.(19), we have after some algebra 



where 



and 



■'^^ s+i(f2^-tJo)+7u 



A{s) = / dtA{t)e''' (2^ 

Jo 



00 

St 



Zis) = / dtZ(t)e 
Jo 



or, 



Zis) = -rY^g^CM ^ , (29) 



where we have used the fact that the amplitudes and the phases of C^(0) are random. Also we 
have A(0) = a(0). Here we have chosen the initial time to = for convenience. 

We make Wigner-Weisskopf approximation to solve for the zeros of A in Eq.(27), where 
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For weak interaction, zeroth approximation is A = if s = 0. As a next approximation let 
s — > in the denominator of the sum in A. In other words under Wigner-Weisskopf approximation 
we calculate the first order shift in the simple pole which is approximately given by 



A (0) - s ^ Lt_o E ^-fO^'" ^^ c = 7^ + ^^^ (31) 

where 7'^ and 5uj are real quantities. Explicit calculation in the usual way yields, 



7^ = / dnp{^)g\n)- ^-^— ^, (32) 

and 



5u = fdnpin)g'in) J^ y| ^ . (33) 

The expressions for the line width 7*^ and the frequency shift 6uj of the system mode thus 
obtained due to the relaxation of the nonequilibrium bath are markedly different from the usual 
expressions of Wigner-Weisskopf theory. This is because of the explicit dependence of 7^^ and 600 
on the 7 in Eqs. (32) and (33) which arises due to the relaxation of the nonequilibrium modes due 
to their coupling with the thermal bath. In the limit 7^0 one recovers the usual decay rate F 
and the level-shift terms, i.e.. 



lim7^ = F . 

7->0 

We thus see that the effect of Wigner-Weisskopf approximation is to replace the more exact 
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equation(19) by the Langevin equation whose solution obtained after appropriate inverse Laplace 
transform of equation(27) is given by 



A{t) = a(0)e-(^"+^^-)*-^(?,C,(0) 



7t„-i(Sl;j-t<j())t 



e ' e 



]_ _ g-(7^-7)tg«(f^M-'^0-5a;)t 



{uo -n^ + Sou) - ^(7^ - 7) 



(34) 



It may be noted that while deriving 7^ and 600, weak coupling and smooth density of states 
p(r2) for the intermediate oscillators have been assumed. By the same token, the correction terms 
in Eq. (34) are small and the result for A{t) as expressed in Eq.(34) could be simplified further 
to the following form, 



A(t)=a(0)e-(^"+^^-)* + x(cuo)a„(0) 



g-7t _ ^-{■y'^+iSuj)t 



(35) 



where x(^o) is given by 



xM 



g{ujo)p{uJo) 
6u - i (7^ - 7) 



The primary evolution of the system mode A thus depends on the secondary relaxation of the 
intermediate oscillators explicitly. 

The expression for 7^ (Eq.(32)) illustrates a dynamical modification of the Wigner-Weisskopf 
decay rate constant since it incorporates the effect of coupling of the nonequilibrium bath to the 
thermal bath through 7. The modification of atomic spontaneous decay rate both in the form of 
enhancement and suppression by appropriate tailoring of vacuum modes of the cavity is now well 
known in cavity QED [7] . Whereas in the cavity QED problems one essentially 77manipulates the 
boundary conditions in various ways, the present modification is effectively dynamical in nature in 
the sense that it carries the effect of relaxation of the nonequilibrium modes on the dissipation of 
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the system mode. It is then also expected that the atomic decay rate might be similarly affected 
in appropriately modified situations. 

Before completing this section we point out that in the present problem of quantum theory 
of dissipation in the quantum optical context we are concerned with the frequency spectrum of 
the radiation field modes. In the context of solid state the frequency density is assumed to be of 
the Debye type with appropriate regularization by cutoff at high frequency. In the cavity QED 
problems adjustment of boundary conditions may lead to different density of states. Although 
there is no generalization of the dependence of friction on the frequency spectrum, in general, one 
encounters a time-retarded friction. 

V. TIME-DEPENDENT SPECTRUM OF A CAVITY MODE WITH GAIN IN 
CONTACT WITH NONEQUILIBRIUM BATH 

An immediate consequence of the nonequilibrium generalization of the fluctuation-dissipation 
relation is the explicit time dependence of the diffusion constant, as evident in Eq.(26). It is there- 
fore expected that this time- dependence may make its presence felt if one analyses the transient 
noise spectrum of the system mode. With this end in view we now calculate the time-dependent 
spectrum of a cavity field mode coupled to a nonequilibrium reservoir which causes the field mode 
to decay at the rate F. Generally we find the spectrum by applying the quantum regression the- 
orem to the Langevin equation for a quantized field mode interacting with a medium described 
by a gain a. In general, the complex gain a(t) is an operator that is saturated by the number 
operator J^(i)A(i) [p. 467 of Ref.l2]. The Langevin equation for our problem is given by, 

A{t) = -(r + i6- a)A{t) + Z{t) , (36) 

where A{t) and Z{t) are slowly varying annihilation and noise operators, respectively. In general, 
Eq.(36) applies to laser- like situations including those with two-level and semiconductor media 



[12]. Here S{= ujo — v) is the detuning of the mode oscillation frequency v from the passive cavity 
resonance frequency cjq and a, the gain coefficient is assumed to be a real number (F > a). Z{t) 
is the noise source operator as given by 

Z{t) = -lY^aACtit) + C^(to)e(-^^''-^'^-^-)(*-*«)]. (37) 

The noise is characterized by the following properties : 

{Z{t))^^ = 0, (38) 

(Zt(t)Z(t'))^^ = [rn + {D{to) - rn}e-2,(i-,,)] ^(^ _ ^/) _ (39) 

It has also to be noted that since we are dealing with a non-stationary situation the standard 
steady state spectrum is not applicable. We therefore take resort to non-steady state spectrum 
or the so-called "physical spectrum" of the cavity mode [8] where the attention is focused on the 
dynamic evolution of the spectrum following an abrupt excitation of a near-resonant cavity mode. 
The main reason for studying the time-dependent spectrum is that the familiar power spectrum of 
the Wiener-Khintchine theorem is not applicable to nonstationary processes. Although in quite a 
number of earlier cases time-dependent spectrum of Page and Lampard were used widely serious 
objections were raised against this spectrum (e.g., it can be negative). Eberly and Wodkiewicz 
have shown that the suitably normalized counting rate of a photodetector can be used to define 
a time-dependent spectrum or physical spectrum. This allows the infiuence of the spectrum 
analyzer (basically a Fabry-Perot interferometer, for example) to be exhibited in the spectrum 
so that the band limit of the measuring device is appropriately incorporated which makes the 
spectrum free from ambiguities and unphysical characteristics of earlier spectrum. It has also 
been emphasized [9] that when the instrument width, W is narrow enough such that W^ <^ F , 
the spectrum appears to be qualitatively similar to Wiener-Khintchine spectrum. This transient 
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spectrum has been used earlier in several occasions in connection with resonance fluorescence 
studies [9], micromaser problem [10] etc. One can define the time-dependent spectrum or the 
physical spectrum as follows; 



S(t, u, W) = 2WRe I dt2e-'^^'-'''^ f ' dre^^"'^^^ {A\t2 + r)A{t2)). (40) 

JO Jo 

Here the symbols have the following meaning : t is the elapsed time after the system and 
the reservoir have been subjected to initial excitation at t = tQ{= 0),W is the full width of the 
transmission peak of the interferometer and A(= uj — u) is the detuning, or frequency offset of 
the Fabry-Perot line center above the frequency of the field uj. It is important to note that the 
time-dependent spectrum is expressed in terms of two integrals in Eq.(40). The first integral is 
over the correlation time and is actually the counterpart of the Wiener-Khintchine spectrum band 
limited by the width of the measuring device, W, while the second one over t2 takes into account 
of the nonstationarity which makes the correlation function t2 dependent. The device width in 
the second integral also sets the limit over the time scale of this nonstationarity. 

Since the time-evolution of the system is governed by Eq.(36), one can make use of the quantum 
regression hypothesis which yields two-time correlation function 

{A\t + T)A{t)) = e-^^-'^"'^^A^{t)A{t)). (41) 

We emphasize that t in Eq.(41) [or t2 in Eq.(40)] refers to the non-stationary time. We therefore 
calculate the explicit time dependence of {A''{t)A{t)) using Einstein's relations (see Appendix A 
for details); 



^{A\t)A{t)) = -2(r - a){A\t)A{t)) + 2rn[l + re-'^*] , where r = ^^ - 1 . (42) 

at D[oo) 

The solution of Eq.(42) is given by 

20 



{A\t)A{t)) 



-2{V-a)t 



(At(O)A(O)) 



+Vn(ujQ) 



;i+r)(r-a) 



V — a (r — a)(r — a — 7) 



7g-2(r-a)t^^ 



-27* 



(r - a - 7) 



(43) 



From Eq.(42) we obtain the steady state condition 



{T - a){A\oo)A{oo)) =Tn. 



(44) 



We let {A^ {00) A{oo)) = N. Eqs.(44) and (43) may then be rewritten, respectively, as 



r — a ,^ 

n= „ A^ 



(45) 



and 



{A\t)A{t)) = N{1 - rke^^''' + rke-^^'), 



(46) 



where a = F — a, /c = 



a— 7 



Combining Eqs.(41) and (46) we obtain the two-time correlation function 



{A^{t2 + T)A{t2)) = e-(^-*'^-")"iV(l - rke-^""'^ + rke'^^'^). 



(47) 



Making use of this relation and performing the integration over r and ^2, we extract the real 
part [Eq.(40)] which yields 
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S{t,A,W) 



2NW 

iy2 + A2 



W- 



^^^ . a + A;r)e-2^* - -kre-'^"^^^^' - (1 + -kr)e- 



Wt 



W+W- + A^ 



|l + /ere ^'^*| + /ere" 



2at 



;i + A;r)2aA 

W^2+A2 ' 



W+t 



sin At 



2afcr - {W^W^ + A2) _^^^ 
■^ W^2 + A2 "^ 



cos At 



(4J 



where, W+ = ^ + {T-a) and 1^_ = f - (F - a) . 
Here we have set the detuning 5 = 0. 
In the long hmit (t -^ oo) the spectrum reaches the steady state value 



S{/\,W) 



2NW 



{f-(r-«)} +A^ 



f)' 



fr 



a 



A2 



{f + (r-a)} +A2 



(49) 



It is interesting to note that at short time the effect of nonequilibrium bath is prominent through 
7 and r. While r(= ^M — 1) includes the effect of preparation of the initial nonequilibrium 
condition by a sudden external excitation at t = which makes the initial diffusion coefficient -D(O) 
to be different from its stationary long time value D{oo), 7 carries the information of relaxation. As 
expected, the steady state spectrum is independent of both 7 and r. This is because at large time 
when the nonequilibrium bath is equilibrated, the system forgets its past and the time-dependence 
of diffusion coefficient is erased and the spectrum becomes the steady state spectrum. 

We now look for the transient characteristics of the spectra of the cavity mode. In Fig 1 
we plot the spectra for different scaled time t after the initial excitation, F being used as a 
scaling parameter. For numerical computation we choose the following scaled parameter set : 
W = 4,a = 0.1,r = 1 and 7 = 0.1. One observes that after the initial excitation the spectra grow 
to a maximum and then the peak height start diminishing and eventually reaches the steady state 
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value. Thus the effect of relaxation of the nonequilibrium bath becomes prominent in the short 
time region. The variation of peak intensity with time for the three different 7 has been shown in 
Fig2. It is apparent that the spectrum reaches the steady state more quickly for larger values of 
7 and also for small 7 the maximum peak is larger than that for larger 7. 

VI. CONCLUSIONS 

We have developed the quantum theory of dissipation of a harmonic oscillator coupled to a 
nonequilibrium bath in terms of a microscopic model. Making use of appropriate separation of 
time scales one can construct an effective Langevin dynamics with memory (where the memory 
functions are not the phenomenological inputs) which is due to the relaxation of the nonequilibrium 
bath modes and identify the relevant noise sources. An essential offshoot is the nonequilibrium 
generalizations of the familiar fluctuation-dissipation and Einstein's relations. It is important to 
note that the Wigner-Weisskopf decay rate constant of the oscillator is dynamically modified. The 
theory is further applied to calculate the time-dependent spectrum of a cavity mode with suitable 
gain. One observes that the nonequilibrium nature of the bath modes makes its presence felt in 
the time-dependence of the diffusion constant leading to the differential behavior of the transient 
spectra from the steady-state ones. Although in the present problem we are mainly concerned 
with the dissipation of energy, it may also be worthwhile to investigate the problem of decoherence 
on similar footing. We hope to address this and related issues in future. 

Acknowledgment : Thanks are due to the Department of Science and Technology (Govt. of 
India) for partial financial support. 
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Appendix A 
Generalized Einstein's relations 

In this section we outline the derivation of the nonequilibrium generalization of Einstein's 
relations Eq.(41). 

The system operators follow the Langevin equation of motion Eq.(35) 



A = -{T + iS-a)A + Z{t) ; (50) 

where the first term within the parenthesis of the right hand side is the drift term and Z{t) is the 
nonequilibrium noise operator [Eq.(37)] . From the identity 

A\t) = A\t - At) + r dt'A^t') (51) 

Jt-At 

we first obtain the correlation function of the system and the noise operator ; 

{A^{t)Z{t)) = {A^{t - At)Z{t)) + f dt\[-{T -id- a)A\t') + Z{t')\Z{t)). (52) 

Jt-At 

Because the operator A'^{t') at time t' is not affected by fluctuation at a later time t, the first 
term on the right hand side is zero. Similarly the correlation {A\t')Z{t)) is zero except at the 
point t' = t but the integration is zero. Thus we have 

{A\t)Z{t)) = r dt\Z\t')Z{t)). (53) 

Jt-At 

Note that we have not assumed the stationary property of the noise. 
Next we determine the equation of motion for the average (A^(t)A(t)); 
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j^{A\t)A{t)) = {A\t)Ait)) + {A\t)A{t)). 



(54) 



From Eq.(50) we have after some algebra 



dt 



{A\t)A{t)) = -2(r - a){A\t)A{t)) + {Z\t)A{t)) + {A\t)Z{t)). 



(55) 



Substituting Eq.(52) and its adjoint in Eq.(55) and performing the integral over t' where we 
use Eq.(39), we obtain the nonequilibrium generalization of Einstein's relation 



dt 



{A^{t)A{t)) = -2(r - a){A^{t)Ait)) + 2Tn{coo) 



1 + 



D(oo) 



1 I e"^'^^*"*'') 



(56) 



where we have used the notation D{oo) = Tn{ujQ). 
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Figure Captions 

1. Time-dependent spectra of the cavity mode witli gain for different dimensionless times 
(a) rt = 0.05, (b) rt = 0.5, (c) Tt = 3.0, (d) Tt = 10.0 with a = 0.1, 5 = 0.0, W = 4.0 and 
r = 1.0 (Scales arbitrary). 



2. Variation of intensity of the time-dependent spectra vs. dimensionless time for different 
dimensioless decay rate constants of the nonequilibrium modes; 
(a) ^ = 0.1, (b) f = 0.2, (c) Y = 0.3 with a, 6, W and r same as Fig. 1 (Scales arbitrary). 
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